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ABSTRACT 

The  distortion  of  flames  in  flows  with  vonical  motion  is  examined  via  asymptotic  analysis  and 
numerical  simulation.  The  model  consists  of  a  constant-density,  onc-siep,  irreversible  Arrhenius 
reaction  between  initially  unmixed  species  occupying  adjacent  half-planes  which  are  then  allowed 
to  mix  and  react  in  the  presence  of  a  vortex.  The  evolution  in  time  of  the  temperature  and  mass- 
fraction  fields  is  followed.  Emphasis  is  placed  on  the  ignition  time  and  location  as  a  function  of 
vortex  Reynolds  number  and  initial  temperature  differences  of  the  reacting  species.  The  study 
brings  out  the  influence  of  the  vortex  on  the  chemical  reaction.  In  all  phases,  good  agreement  is 
observed  between  asymptotic  analysis  and  the  full  numerical  solution  of  the  model  equations. 


Ihc  .second  and  third  authors  were  supjxined  hy  the  Nalii)n.il  .'\eroiiauucs  and  Space  Admirnstraliun  under  N.^S.\  Comracl  No.  N.ASl- 
IS605  while  in  residence  at  the  Insiiiule  lor  Computer  Applicalions  in  Science  and  lingmeeniii!,  NASA  Langley  Research  Cenicr. 
Hampton,  VA  23f)()5. 
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1.  INTRODUCTION 


In  a  seminal  paper.  Marble  (1985)  developed  a  simple  model  problem  which  describes  the  distortion 
of  a  diffusion  flame  by  the  presence  of  a  viscous  vortex  flow,  as  might  be  found  in  the  large  scale  sffuciures 
of  turbulent  diffusion  flames.  Assuming  fast  chemical  kinetics  and  the  constant  density  approximation,  a 
flame  sheet  exists  at  time  zero  separating  a  region  of  fuel  in  the  upper  half  plane  and  a  region  of  oxidizer  in 
the  lower  half  plane.  As  time  increases,  the  flame  sheet  is  allowed  to  evolve  by  diffusion  and  convection 
and  is  wound  up  by  the  vortex  flow  field.  The  analysis  was  based  on  the  assumption  that  the  flame  elements 
are  locally  strained  and  sheared,  but  that  the  shearing  motion  was  unimportant  and  hence  could  be 
neglected.  Theoretical  analysis  shows  that  (i)  the  vortex  field  increa.ses  the  reactant  consumption,  (ii)  the 
reaction  rate  is  independent  of  time,  (iii)  a  similarity  rule  for  the  core  radius  growth  exists,  and  (iv)  a  simi¬ 
larity  rule  exists  for  the  reactant  consumption  rate. 

The  analytical  problem  posed  by  Marble  (1985)  (henceforth  called  the  Marble  problem)  has  been 
solved  numerically  by  Laverdant  and  Candel  (1988)  and  Rehm,  Baum,  Lozier  and  Aronson  (1989),  but 
without  the  strain-shear  approximation.  Thus  the  complete  flow  field  was  considered.  Laverdant  and  Candel 
(1988)  solved  the  Marble  problem  by  an  implicit  finite  difference  scheme.  They  verified  numerically 
Marble’s  similarity  rule  for  the  core  radius  growth  and  the  result  that  the  reaction  rate  is  independent  of 
time,  and  also  verified  the  correemess  of  the  shear  approximation.  In  addition,  the  effect  of  varying  the 
stoichiometric  (or  equivalence)  parameter  was  also  studied.  They  showed  that  die  core  is  mostly  surrounded 
by  the  oxidizer  when  the  stoichiometric  parameter  is  less  than  unity  (fuel  lean),  and  is  mosUy  surrounded  by 
the  fuel  when  it  is  greater  than  unity  (fuel  rich).  Rehm,  et  al.  (1989)  solved  the  Marble  problem  by  first 
transforming  the  system  into  a  Lagrangian  coordinate  system.  Then  assuming  a  similarity  solution,  the 
resulting  problem  was  solved  analytically  using  Fourier  series.  The  Fourier  amplitudes  were  determined  by 
either  a  two-point  boundary-value  solver,  valid  for  all  values  of  the  independent  parameters,  or  in  the 
asymptotic  limit  of  large  Schmidt  number  (as  might  be  more  appropriate  for  liquids  than  for  gases).  The 
results  of  this  study  confirmed  the  behaviour  of  the  numerical  solution  of  Laverdant  and  Candel  (1988),  as 
well  as  the  dependency  of  the  reactant  consumption  rate  found  by  Marble  (1986).  Norton  (1983)  also 
solved  the  Marble  problem  with  the  strain-shear  approximation,  but  using  finite  rate  chemistry  instead  of  the 
flame  sheet  approximation,  confirming  the  results  of  Marble. 

Many  extensions  of  the  Marble  problem  exist,  each  utilizing  the  same  analytical  approach  as  that  of 
Marble  (1985)  and  assuming  the  fast  chemistry  limit,  .some  of  which  are  briefly  reviewed  here.  Karagozian 
and  Marble  (1986)  considered  the  time  dependent  interaction  of  a  diffusion  flame  with  an  axially  strained 
vortex  in  the  third  dimension.  It  was  .shown  that  the  main  conclusions  for  the  three  dimensional  vortex  field 
remains  formally  identical  to  die  two  dimensional  vortex  held  (the  Marble  problem).  In  addition,  density 
effects  were  shown  to  be  small.  Karagozian  and  Manda  (1986)  considered  the  distortion  of  a  two  dimen¬ 
sional  fuel  strip  in  the  presence  of  a  pair  of  counter-rotating  vonicies.  The  fuel  strip  was  assumed  either 
infinite  or  semi-infmitc  in  extent,  .separated  from  the  oxidizer  by  two  flame  sheets.  On  each  flame  sheet 
resides  a  vortex,  which  are  allowed  to  rotate  in  opposite  directions.  The  .semi-inlinite  fuel  stnp  is  particu¬ 
larly  relevant  to  the  vortical  flame  structure  formed  at  a  circular  ontice  or  nozzle.  Cetegen  and  Sirignano 
considered  the  interaction  of  a  diffusion  flame  with  a  single  vortex  (1987)  and  the  interaction  with  an 
infinite  row  of  two  dimensional  vortices,  which  is  representative  of  the  temporal  growth  of  a  mixing  layer 
before  pairing  has  occurred  (1988).  In  both  cases,  die  emphasis  was  on  the  construction  of  the  concentration 
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profiles  in  the  vortical  structure,  tis  well  as,  the  construction  ol'  the  probability  density  functions  (pdfs).  A 
review  of  these  two  cases  can  be  found  in  Sirignano  (1989).  Peters  and  Williams  (1988)  developed  an 
analogous  problem  for  a  premixed  Ihune,  relevant  for  describing  premixed  turbulent  combustion  in  large 
scale  suuctures.  The  emphasis  of  this  study  was  in  describing  the  growth  of  the  reacted  core,  flame  extinc¬ 
tion  by  suetch  for  Lewis  numbers  greater  than  unity,  and  the  effect  of  heat  release. 

In  iuiothcr  seminal  paper,  Linan  and  Crespo  (1976)  examined  the  structure  of  a  diffusion  flame  in  the 
unsteady  mixing  of  tw'o  half  spaces  of  fuel  tuid  oxidizer.  By  using  a  combination  of  large  activation  energy 
asymptotics  and  numerics,  tliey  analyzed  the  continuous  evolution  from  nearly  frozen  flow  to  near  equili¬ 
brium  flow  and  showed  that  three  laminar  regimes  exist  within  the  flow;  the  ignition,  deflagration,  and 
diffusion  flame  regimes.  The  ignition  regime  is  a  region  where  the  combustible  gases  mix  until,  at  some 
finite  time,  a  thermal  explosion  occurs  at  a  W'ell  delined  location  and  the  gas  is  ignited.  The  second  regime 
is  the  deflagration  region.  After  ignition,  a  pair  of  deflagration  waves  (or  "premixed  flamcleis")  emerge 
according  to  classical  thermal  explosion  theory.  These  waves  ari.se  because  the  mixture  is  pot  stoichiometric 
in  the  premixed  region.  One  of  tlie  reactants  is  consumed  locally,  leaving  behind  an  excess  of  the  other 
reactant.  Thus,  one  of  the  (lame lets  is  luel-rich  and  the  other  is  fuel-lean.  There  is  excess  fuel  behind  the 
fuel -rich  llamclet  and  excess  oxidiz.er  behind  the  fuel-lean  llamelet.  Concentration  gradients  behind  the 
llamelcLs  drive  the  unburnt  fuel  and  oxidizer  towards  the  diffusion  flame  where  they  are  consumed  locally. 
These  premixed  llamelets  are  quite  weak  in  that  the  temperature  ri.se  associated  with  them  is  small,  and  they 
exist  only  until  all  of  the  delicient  reactant  is  consumed.  Just  beyond  the  deflagration  w-aves,  a  diffusion 
flame  regime  exists  where  the  mixing  process  is  governed  by  diffusion  in  the  direction  normal  to  the  flame. 
As  time  increases,  the  diffusion  flame  approaches  a  flame  sheet.  We  note  that  the  existence  of  a  well 
defined  ignition  point  and  the  premixed  llameleis  depends  critically  on  tlie  relative  magnitudes  of  the  two 
initial  temperatures  to  that  of  tlie  adiabatic  diffusion  flame  temperature.  That  is,  if  the  adiabatic  diffusion 
flame  tenifx'rature  is  greater  than  either  inititil  temperature,  a  well  defined  ignition  point  always  occurs,  fol¬ 
lowed  by  the  premixed  llamelets.  On  the  other  hand,  if  the  adiabatic  diffusion  flame  temperature  is  between 
the  two  initial  temperatures,  there  is  no  w-ell  defined  ignition  point,  and  a  single  premised  flame  merges 
smoothly  into  the  diffusion  fliuiie.  The  analytical  results  of  Linan  and  Cresjx)  were  extended  by  Jackson  and 
Hussaini  (1988)  and  Gro.sch  and  Jackson  (1991)  in  their  studies  of  steady  supersonic  combustion  in 
compressible  laminar  mixing  layers.  In  particular,  the  three  flame  regimes  also  exist  at  high  speeds  provided 
the  adiabatic  diffusion  flame  temperature  is  greater  than  cither  freestream  temperature.  Finally,  we  note  that 
this  configuration  of  two  llamelets  and  a  diffusion  flame  (collectively  called  a  u-ibrachial  flame)  has  also 
been  observed  m  the  propagation  of  a  premixed  flame  into  a  nonuniform  mixture  (e.g.,  Buckmaster  and 
Maialon,  1988). 

Our  goal  is  to  combine  the  imxiels  of  .Marble  (198.5)  and  Linan  and  CTc.sjX)  (1976).  That  is,  we  pro¬ 
pose  to  investigate  the  continuous  evolution  of  initially  unmixed  species  occupying  adjacent  half  planes  to 
near  equilibrium  flow  in  the  presence  of  a  viscous  vortex  without  invoking  the  fast  chemistry  limit. 
Emphasis  is  placed  on  the  ignition  time  and  location  as  a  function  of  vortex  Reynolds  number  and  initial 
lempeialure  difterences  of  the  reacting  species.  In  the  next  section  we  slate  the  problem.  .Section  3  presents 
an  analysis  ol  the  ignition  regime  using  a  combination  of  large  aeiivatiun  energy  fund  hence  Zeldovich 
number)  asymptotics  and  numerics.  Scvtion  4  presenus  a  simple  model  problem  for  the  diffusion  flame 
regime  describing  the  viscous  core  region  about  the  flame  sheet,  and  identifies  certain  physical  effects  by 
employing  asynipioiic  expansions  and  numerical  solutions.  The  full  equations  are  tlien  .solved  numerically 
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in  section  5.  Comparisons  between  asymptotics  and  numerics  are  given  through  out  the  study.  Finally,  our 
conclusions  are  presented  in  section  6. 

2.  EQUATION  FORMULATION 

In  this  section  the  problem  for  the  time  evolution  of  initially  unmixed  species  occupying  adjacent 
half-planes  which  are  then  allowed  to  mix  and  react  in  the  presence  of  a  vortex  is  stated.  The  nondimcn- 
sional  equations,  boundary  conditions  and  initial  conditions  governing  this  field  in  cylindrical  coordinates 
are 

U=0.  F  =  ^^(1  -  cxp(-r2/4/ 5c)),  (2.1a) 

r 

7',  + -^7'e  =  V‘T  +  BQ.  (2.1b) 

r 

=  7  =  1,2.  (2,1c) 

Q  =  DriF,F2e'^  ^  .  (2. Id) 

7’  =  =  1,  ^2  =  0  at  /  =  0,  r  >  0,  0  <  0  <  7t,  and  /  >  0,  r  — »  oo,  0  <  6  <  n,  (2.2a) 

^2- 

T  =  P/ ,  F 1  =  0,  f  2  =  —  =  <t>  tit  t  -  0,  r  >  0,  K  <  i)  <  2k  , 

and  t  >  0,  r  — >  oo,  71  <  0  <  2  71,  (2.2b) 

where  (U,V)  are  the  radial  and  angular  velocity  com|x)nents,  respccuvely;  is  the  two-di-  .  nsional  Lapla- 
cian  operator  in  cylindrical  coordinates;  T  is  the  temperature;  and  F,  and  F2  the  mass  fr  ciions.  The  chemi¬ 
cal  model  is  a  one-step  irreversible  Arrhenius  reaction,  and  the  vortex  model  is  an  Oscen  vortex  with  circu¬ 
lation  r.  The  nondimensional  parameters  appearing  above  are  the  vortex  Reynolds  number  R ,  dclined  as 
r  /  2tcv  where  v  is  the  kinematic  viscosity;  the  Schmidt  number  5c  =  v/D  a*^^ aimed  equal  for  both  species 
and  D  is  the  species  diffusion  coefficient;  the  Zeldovich  number  Ze  =  E  /''V'  T^)  with  F  the  dimensional 
activation  energy  and  R"  the  universal  gas  constant:  the  Damkohicr  nun. her  Da  defined  as  the  ratio  ol  the 
characteristic  diffusion  time  scale  to  the  characteristic  chemical  time  ,cale;  p  the  heat  release  per  unit  mass 
of  Fj;  the  equivalence  ratio  4)  defined  as  the  ratio  of  the  mass  fmetion  Fj  to  the  mass  fraction  F 2  divided 
by  the  ratio  of  their  molecular  weighes  times  their  stoichiometric  coefficienLs;  and  finally  P^  the  ratio  of  the 
initial  temperatures.  If  ())  =  1,  the  mixture  is  said  to  be  stoi  liiometric,  if  i))  >  1  it  is  F,  rich,  and  if  (j)  <  1,  it 
is  F|  lean.  Also,  if  B.  is  less  than  one,  F2  is  relatively  cold  compiled  to  F,,  and  if  P/  is  greater  than  one  it 
is  relatively  hot.  The  temperature  and  mass  fractions  were  nondimensionali/ed  by  the  values  and  F  < 
respectively.  The  velocities  were  nondimcnsionali/ed  by  f/„.  some  characteristic  six'cd.  LengtJis  and  times 
arc  referred  to  the  relevant  diffusion  chciacteristic  scales  --  X/  and  Ij/U^  of  die  How, 
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respeciively,  where  p  is  the  density  assumed  constant,  X  is  the  thermal  conductivity,  and  Cp  is  the  specific 
heat  at  constant  pressure.  Finally,  tlie  assumption  of  unit  Lewis  number  is  made,  so  that  Sc  =  Pr,  where  Pr 
is  the  Prandtl  number. 


3.  IGNITION  REGIME 


At  time  f  =  0,  the  reaction  rate  is  exactly  zero  owing  to  the  product  F1F2  =  0.  For  t  >  0,  the  fuel  and 
oxidizer  begin  to  mix  by  diffusion,  as  well  as  by  convection  due  to  the  presence  of  the  vortex,  and  the  reac¬ 
tion  rate  is  no  longer  zero.  For  small  time,  it  can  reasonably  be  assumed  that  the  effect  of  the  reaction  is 
small,  and  in  the  absence  of  reaction  term  the  system  (2.1)  reduces  to 

7;  +  ^  (1  -  e.xp(-r2/4i  Sc)]7'6  -  V^T,  (3.1a) 

r 


-^[1 -cxp(-r‘/4r5c)JF,e=  7  =  1,2.  (3.1b) 

The  above  system  must  be  solved  numerically  subject  to  the  conditions  (2.2).  and  its  solutions  are  referred 
to  as  the  inert  solutions,  denoted  by 

7)ay)(r.0,O,  Fjj  j  =  \,2.  (3.2) 

As  time  increases,  more  of  the  combustible  mixes  until,  at  some  finite  time,  a  thermal  explosion  occurs 
characterized  by  signilicant  departure  from  the  inert. 

To  analyze  the  ignition  process,  we  determine  the  effect  of  the  growing  reaction  rate  by  expanding 
about  the  inert  solution  as 

T  =  T,  +  Ze-'  Tu  F,  =  F,j  +  7  -  1,2,  (3.3) 

and  take  the  asymptotic  limit  Ze  Substituting  into  (2.1),  we  sec  that  the  left  hand  sides  are  0(Ze“‘) 

while  the  right  hand  sides  arc  0{Da  e  '  ’),  where  7’.  =  tmxQ'i).  In  order  to  balance  these  terms.  Da 

must  therefore  have  the  form 


„  1  Zf/T. 

'  pr' 


Substituting  into  (2.1b)  yields  the  following  equation  for  T’j 


7'i.,  +  ^  11  -  cxp(-r2/4i  .'ic)17-,,e  = 


/  1  -f  {I I  —  7 •  ) 

17(1)  +  Zc^'F,)” 


(3.4) 


+  lFi_,  +  Ze-'F^,,  I  IF.,,  +  Zc-‘F.  ,  1  exp 


(3.5) 
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which  must  be  solved  subject  to  the  initial  and  boundary  conditions 

7'i  =  0  at  t  -  0,  r  >  0,  0  <  Q  <  2n,  and  t  >  0,  r  ->  £»,  0  <  0  <  2 7t.  (3.6) 

Once  Ti  is  known,  the  mass  fraction  perturbations  can  be  found  directly  from  (2.1c).  Note  that  the  right 

hand  side  of  (3.5)  is  exponentially  small  except  where  I'/)  -  T.  j  =  0(Ze“').  There  are  now  two  cases  to 

consider,  depending  upon  the  magnitude  of  the  parameter  1 1  -  Pj-  | .  In  this  section  we  only  consider  the  case 

of  nearly  equal  initial  temperatures  (P^  =  1),  and  will  postpone  the  analysis  for  0(1)  initial  temperature 
differences  as  a  topic  for  future  research. 

As  mentioned  above,  only  the  important  case  of  ignition  for  nearly  equal  initial  temperatures  of  the 
fuel  and  oxidizer  is  considered.  To  this  end,  we  set 


p,  =  l+Ze-'Pr.  (3.7) 

with  Py  fixed  and  0(1).  With  this  choice,  die  inert  temperature  7)  is  given  by 

7/  =  1  +  Ze ,  (3.8) 

where  7)  j  satisfies  (3.1a)  subject  to  the  initial  and  boundary  conditions 

7/1  =  0  at  t  =  0,  r  >  0,  0  <  6  <  K,  and  t  >  0,  r  0  <  0  <  n,  (3.9a) 

^/.i  =  ^r  r  =  0,  r  >  0,  rt  <  0  <  2k,  and  t  >  0,  r  ->  «,  k  <  0  <  2k.  (3.9b) 


Substituting  into  (3.5),  together  with  the  asymptotic  choice  7.  =  1,  yields  the  following  equation  for  the  dis¬ 
turbance  temperature 

T,.,  +  ^[1  -exp(-r^/4rSc)]7-,  +  (3.l0) 


subject  to  the  conditions  (3.6).  This  equation  was  solved  numerically  for  a  range  of  7?  and  Py  with 
Sc  -  0  =  1.  Implementation  of  boundary  conditions  arc  facilitated  if  the  system  is  recast  in  Cartesian  coor¬ 
dinates.  The  solution  technique  is  a  2nd-order  finite  difference  .scheme  on  a  nonuniform  mesh.  To  resolve 
the  structure  in  the  core  region  of  the  field,  a  coordinate  stretching  of  the  form 


yp  = 


C2-  1 


(3.11) 


is  used,  where  e  [ -yVuu,  >s  the  physical  coordinate,  >■,.  e  1-1,1]  is  the  computational  coordinate, 
and  Cl  and  C2  are  adjustable  constants.  The  smaller  the  quantity  (C2  -  1).  the  su-onger  the  stretching.  In 
the  y  direction,  Ci  is  either  4  or  2,  and  C2  is  2.  A  similar  sU'etching  is  used  in  the  x  direction,  w'ith  Ci  =  4 
and  C2  =  2.  To  avoid  the  singularity  at  the  origin,  no  mesh  points  are  placed  there.  The  outer  boundaries 
are  set  at  .50  or  2(X)  for  Xn,,,,  and  20  or  50  for  yma*-  Grid  resolution  studies  which  at  least  doubled  the  com¬ 
putational  mesh  were  carried  out  to  ensure  that  structures  were  well  re.soKed.  The  resolutions  required 
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ranged  from  a  64^  mesh  lo  a  256^  mesh  for  large  vortex  Reynolds  number  (R  =  5000).  The  time-stepping 
scheme  is  a  four-stage  Rungc-Kutla  which  is  formally  2nd-order  but  has  an  extended  stability  region  making 
it  accurate  and  robust  for  moderately  stiff  problems.  All  runs  were  performed  interactively  on  a  Cray  YMP. 

As  t  increases,  the  solution  for  1\  becomes  unbounded  at  some  finite  time  (r.^)  and  location 
This  characterizes  the  ignition  regime.  The  special  case  R  =  0,  in  which  two  initially  unmixed  species  are 
allowed  to  diffuse  without  the  mixing  generated  by  the  vortex,  corresponds  to  the  results  of  Linan  and 
Crespo  (1976).  To  verify  the  solution  technique  described  above,  we  give  a  comparison  of  our  results  to 
that  of  Linan  and  Crespo  in  Table  1.  For  this  case  ignition  takes  place  along  a  line  parallel  to  the  ;t-axis 
which  is  located  at  y  =  0  for  (iy  =  0  and  resides  in  the  hotter  region  for  0. 

The  effects  of  the  vortex  can  be  seen  by  examining  Figures  1  and  2  and  Table  II.  Figure  1  is  a  plot 
of  the  ignition  times  versus  the  vortex  Reynolds  number  R  for  three  values  of  with  Sc  =  (}>  =  1. 
Table  II  contains  selected  ignition  times  as  a  function  of  R  for  the  three  values  of  of  Figure  1.  The 
effect  of  increasing  py  from  zero  (i.c.,  py  >  1)  is  to  enhance  ignition,  while  decreasing  py  from  zero  (i.e., 
Py  <  1)  has  the  opposite  effect.  One  can  see  from  Figure  1  and  Table  II  that  the  vortex  Reynolds  number 
has  little  effect  on  the  overall  ignition  time.  However,  ignition  now  occurs  at  a  point  rather  than  along  an 
entire  line.  This  effect  is  shown  in  Figure  2  where  the  ignition  locations  versus  .r,j  are  plotted  for  vari¬ 
ous  values  of  R  and  P/  =  -2.  Note  that  as  R  increa.sc.s,  the  ignition  location  spirals  clockwise  towards  the 
viscous  core  center.  Ignition  is  seen  to  occur  within  the  core  for  vortex  Reynolds  numbers  R  >  70.  This  is 
in  contrast  to  the  case  of  equal  initial  temperatures  (Py  =  0),  where  ignition  always  occurs  at  the  origin  for 
any  value  of  R .  The  analogous  plot  for  py  =  2  has  essentially  the  .same  characteristics  as  Figure  2.  In  sum¬ 
mary,  ignition  occurs  in  the  region  of  the  initially  hotter  reactant  for  small  R ,  and  spirals  clockwise  towards 
the  viscous  core  center  as  R  increases. 

Figure  2  indicates  that  the  ignition  location  is  at  its  maximum  from  the  center  of  the  vortex  at  approx¬ 
imately  R  =  28.  Surface  ploLs  of  T\  luid  F\  i  shown  in  Figures  3a,b,  respectively,  display  the  skewed  struc¬ 
ture  of  the  reaction  core  at  these  conditions.  The  development  of  the  hot  spot  for  Py  =  -  2  as  a  function  of 
increasing  R  is  shown  in  the  contour  plots  of  Y'l  and  F  u  in  Figures  4  and  5,  respectively.  The  structure  of 
the  reaction  core  develops  from  an  asymmetrically  skewed  center  (Figures  4a-c)  to  an  axisymmetric  center 
(Figure  4d),  However,  the  local  structure  of  the  hot  spot  within  the  reaction  core  is  axisymmeuic  for  any 
value  of  the  vortex  Reynolds  number  greater  than  zero.  This  fact  is  more  cletu’ly  .seen  in  Figure  6  where  we 
plot  /  I  versus  >  at  the  ignition  location  x,g  and  time  t,g  corresponding  to  the  same  conditions  given  in  Fig¬ 
ure  4.  From  this  ligurc  it  is  clear  that  once  the  hot  spot  develops  within  the  reaction  core,  it  is  axisym¬ 
metric. 

For  completeness,  a  conu-ast  of  the  effects  of  a  potential  vortex  distribution  (I'  =  R  Sc  /  r)  versus  the 
Oseen  vortex,  given  in  (2,1a),  on  tiie  ignition  locations  arc  made.  Figure  7  is  a  plot  of  y\g  versus  for  the 
potential  vortex  for  die  .same  conditions  as  in  Figure  2.  Note  that  the  collapse  of  the  spiral  into  the  center 
occurs  for  R  >  7,  as  opposed  to  the  value  R  >  70  for  the  Oseen  vortex. 

The  numerical  solutions  presented  above  suggests  that  near  thermal  runaway,  a  hot  .spot  develops  at 
)  within  the  reaction  core  and  that  the  structure  of  this  hot  spot  is  axisymmetric  (sec  Figure  6).  Thus, 
convection  is  not  important  and  the  local  structure  of  die  hot  spot  is  diffusion  controlled.  In  this  case,  Dold 
(1985.  1989)  has  shown  that  the  proper  .scalings  characterizing  the  local  structure  of  the  hot  spot  tue  given 
by 
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X  =  rig  cos0,g  +  H(T)i  COS0,  y  =  r^g  sin0,g  +  |i(T)i  sin0,  x  =  r,g  -  r 


(3.12) 


where 


la-(x)  T 1  a  -  In  (i2,^  X)]  =  x  4,  (3. 13) 

as  X  — >  0,  with  Q,^  =  j  F2,/  ^  ‘  ^  evaluated  at  {r  ,Q,t)  =  (r,^  ,0,j(  Jig),  and  a  is  a  constant  determined  by 
matching  with  the  initial  conditions  and  is  a  function  of  the  parameters  (R  ,Sc  With  these  scalings,  and 
ignoring  all  dependence  on  0,  the  asymptotic  form  of  7  i  in  the  hot  spot  is  given  by 


T,  = 


X  N  1  2 ,  ,,  ln£  5  +  n  s"IA 

-ln(Q,  X)  -  ln(l  +  .v^/4)  -  — —  - — 

421+  5^/4 


1 

"I 


1  +  n 


5^/4 

1  +^2/4 


ln(l  +  s^l  A) 


+  0{ 


1:2 


(3.14) 


with  n  =  0  if  0  or  /5  =  1  if  r,g  =  0.  Again  note  that  the  local  suuciurc  of  the  hot  spot,  once  it 
develops,  is  independent  of  the  influence  of  the  vortex  and  so  is  diffusion  controlled.  Wc  are  currently 
investigating  the  future  time  development  of  the  hot  spot,  and  the  subsequent  flame  development. 


4.  DIFFUSION  FLAME  REGIME 


Typically,  a  diffusion  tlame  is  characterized  by  a  chemical  reaction  time  that  is  much  smaller  than  a 
characteristic  diffusion  time.  Chemical  reactions  then  occur  in  a  narrow  zone  between  the  fuel  and  the  oxi¬ 
dizer,  where  the  concentrations  of  both  reactants  are  very  small.  Mathematically,  the  assumption  of  very  fast 
chemical  reaction  rates  leads  to  the  limit  of  infinite  Damkohler  number  which  reduces  the  diffusion  flame  to 
a  flame  sheet  (i.e.,  local  chemical  equilibrium).  This  assumption  significantly  reduces  the  complexity  of  the 
problem  since  it  eliminates  the  analysis  associated  with  the  chemical  kinetics.  For  many  flows,  the  assump¬ 
tion  of  local  chemical  equilibrium  adequately  predicts  the  location  and  the  shape  of  the  diffusion  flame 
(Buckmaster  and  Ludford,  1982;  Williams,  1985).  Since  the  flame  sheet  model  is  amenable  to  analysis,  we 
give  its  structure  below.  For  finite  values  of  the  Damkohler  number,  equations  (2.1)  must  be  solved  numeri¬ 
cally  and  this  is  done  in  the  next  section. 

We  begin  the  analysis  of  the  diffusion  flame  regime  by  defining  the  following  conserved  variables 
(see,  e.g.;  Williams,  1985); 


7+pF,  =  P;  +(l-p.^+P)Z. 


(4.1) 


T  +  PF2-P/-  +  P4*'’^(I“Pr“P4>')2^' 


(4.2) 


F|  -  F2  +  ' 

1  +  (J)*’ 


(4.3) 
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where  Z  is  the  mixture  mass  fraction  and  satisfies  the  convection-diffusion  equation 


Z,  [1  -  exp(-/-2/4r  5c)]Ze  =  V^Z,  (4.4) 

sdbject  to  the  initial  and  boundary  conditions 

Z  =  1  at  /  =  0,  r  >  0,  0  <  0  <  K,  and  /  >  0,  r  — >  «>,  0  <  0  <  7i,  (4.5a) 

Z=0  at  t  =  0,  r  >  0,  n  <  Q  <  2n,  and  t  >  0,  r  — >  oo,  71  <  0  <  27t.  (4.5b) 

In  the  limit  of  infinite  Damkohler  number  the  flame  sheet  solution  is  given  by 

=  1  -  (1 -Ki)-')(1 -Z),  F2-O,  (4.6a) 

T  =  p;  +  pr'  +  (1  -  Pr  -  P<t)”')2.  (4.6b) 

valid  for  Z  >  Zy  ,  and 

f,=0.  f2  =  (l>''-(l+(t>'')Z.  (4.7a) 

7=pr  +  (l-p/+P)Z.  (4.7b) 


valid  for  Z  <  Zf .  Here,  Zf  defines  the  location  of  the  flame  sheet  where  both  the  reactants  vanish,  given  by 
the  implicit  relation 


y-f  - 


1  + 


(4.8a) 


and  T  takes  the  adiabatic  flame  value 


1  +  P  y  <l>  +  P 

I  +  (j) 


(4.8b) 


Note  that  tlie  flame  location  is  independent  of  P/  and  p.  Once  Z  =  Z(r,0, y)  is  known,  then  the  other  vari¬ 
ables  (/  .Fi.F;)  can  be  found  from  (4.6)-(4.7). 


As  discussed  m  the  Introduction,  equation  (4.4)  for  Z  lias  been  solved  numerically  by  Laverdant  and 
Candel  (1988)  and  Rehm,  Baum,  Lo/ier  and  Aronson  (1989).  Our  goal  is  not  tlie  numerical  solution  of 
(4.4).  Rather,  we  propose  in  ilie  lollowing  subsections  a  simple  model  problem  describing  the  viscous  core 
center  in  a  region  about  flic  flame  stieei,  wiili  the  goal  of  identifying  certain  physical  effects  by  employing 
asymptotic  expansions. 
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4.1.  LOCAL  ANALYSIS  ABOUT  FLAME  SHEET 


In  this  section  we  formulate  a  model  problem  describing  the  effect  of  the  vortex  in  a  neighborhood 
about  the  flame  sheet.  Assume  at  time  /  =  0  a  flame  sheet  exists  at  y/ .  In  the  absence  of  a  vortex,  the  mix¬ 
ture  mass  fraction  evolves  according  to  (4.4)  with  /?  =  0,  and  the  analytical  solution  can  be  written  in  terms 
of  an  error  function.  At  any  instantaneous  time  ,  the  solution  across  the  flame  sheet  is  easily  seen  to  be  a 
linear  shear  profile  (Figure  8).  given  by 


Z(r,Q,i„)=  +  Z/(/„  )(y  -y/  (to )),  (4.9) 

which  is  only  the  first-order  Taylor  series  expansion  of  the  mixture  mass  fraction  profile  Z  about  the  flame 
sheet  location  y/ ,  found  imphcitly  by  (4.8a).  Here,  the  prime  denotes  differentiation  with  respect  to  y .  The 
model  problem  assumes  that  the  time  scale  associated  with  the  vortex  is  faster  than  that  of  diffusion,  so  that 
on  this  time  scale  the  linear  profile  (4.9)  can  be  considered  steady.  In  addition,  we  assume  that  the  flame 
sheet  resides  close  to  the  origin.  Thus  the  following  scales  can  be  introduced 


r  =  +  £( ,  y  =  -^ty , 


r  =%*£/=, 


=  !-(-%/£$,  Z  =  Y  +  \'£  Z , 


(4.10) 


with  £  <<  1,  and  the  vortex  is  turned  on  at  time  t  =  0.  With  these  scalings  the  linear  shear  profile,  in  terms 
of  Z,  becomes 


Z(r,0,O)  =  -  j  +  Zy'(t<,)[y  -  yy((^)j  s  k^^  +  kiy.  (4.11) 

Note  that  k^  -  -$/4  -  Zf{t^)yf(t^)  measures  the  instantaneous  deviation  of  Z  from  its  sioichiometric  value 
Zy^  =  1  / 2  at  r  =  0,  and  k2  =  measures  the  instantaneous  gradient  of  the  mixture  ma.ss  fraction  across 

the  flame  sheet.  The  convection-diffusion  equation  (4.4)  in  terms  of  the  perturbation  mixture  mass  fraction 
Z  becomes 


Z-  +  ^^[1  -exp(-r^/4/5c)!Ze=  V‘Z,  (4.12) 

r 

where  is  the  Laplacian  operator  in  the  transformed  plane.  The  perturbation  equation  (4.12)  is  to  be 
solved  subject  to  the  initial  and  boundary  condition  given  by  the  linear  shear  profile  (4.11). 

To  examine  the  time  evolution  of  the  flame  sheet  about  the  origin  when  the  vortex  is  present,  we 

define 


Z(r,0,r )  =  I  +  I2I  r  sinB  +  Z(r,0,  f )] ,  (4.13) 

where  Z  is  the  disturbed  mass  fraction  due  to  the  presence  of  die  vortex.  Upon  substituting  (4.13)  into 
(4.12),  we  find  that  Z  satisfies  the  convection-diffusion  equation 
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Z,-+  -exp(-r^/4r5c)](rcosG-hZe)  =  ^^Z,  (4.14) 

f 

subject  to  homogeneous  initial  anti  boundary  conditions.  Note  that  (4.14)  is  non-homogeneous,  and  hence 
does  not  have  a  similarity  solution.  Here,  we  introduce  the  variable 


S' 


(4.15) 


which  is  the  similarity  variable  in  the  absence  o!  the  vortc.x,  so  that  Z  =  Z(s,0,t).  In  the  U’ansfomied  plane 
Z  satisfies 


-  -  •i  -  R  Sr  ,  —  -  -  T  ' 

t  Z-  -  [1  -  exp(-A^'4.Sc)|(5  COS0  -e  Z^)  =  V^Z.  (4.16) 

2  s'- 

This  equation  is  separable  in  space  and  time  and  hence  has  the  simple  analytical  solution 

Z  =  '^t  [Z^  (s)cos0  +  Z^  (.v)h:in0],  (4.17) 

where 

\(Z,  -sZ;)+  ^11  -cxp(-s2/45c)](Z„  +s)^z;’  +  -z;-  \z,  ,  (4.18) 

2  5  5 

4(Z.  -iZj)- -^^(1 -cxp(-s2/45c)iz,  =t"+ -t'-  .  (4.19) 

subject  to  the  boundary  conditions 

2,  (0)  =  Z„  (0)  =  Z,  («)  =  Z„  (-)  =  0.  (4.20) 


Here,  primes  denote  ordinary  ddlercntiation  with  respect  to  the  similarity  variable  s,  Z^  is  the  amplitude  of 
the  even  mode,  and  Z^  is  the  amplitude  of  the  odd  mode.  Note  that  the  above  solution  is  valid  for  ail  times 
and  vortex  Reynolds  numbers. 

The  system  (4.18-4.19)  subject  to  (4.20)  was  solved  numerically  using  a  2nd-order  finite  difference 
scheme  with  an  appropriate  sU'ctching  in  s.  The  solution  to  this  system  for  a  range  of  vortex  Reynolds 
numbers  is  given  in  Figure  9.  As  the  vortex  Reynolds  number  increases,  the  number  of  oscillations  in  each 
of  the  components  Z,  and  Z„  also  increases  while  mo'  ing  away  from  the  origin,  thus  establishing  a  core 
region.  The  composite  solution  (4,13)  is  plotted  in  Figure  10  for  t  =  1,  0  =  ±Tt/2,  k2=  1,  itj  =  0,  and 
5i.'  =  1.  From  this  figure  the  extent  of  the  core  region  is  seen  to  grow  as  R  increa.scs,  and  this  core,  once 

pr 

established,  grows  like  (sec  (4.17))  as  noted  previously  by  Marble  (1985). 
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4.2.  LARGE  TIME  BEHAVIOUR 

To  investigate  the  local  structure  of  the  model  problem  described  in  the  previous  section  for  large 
times,  let  t  ->  <»,  so  that  s  ->  0  for  fixed  r.  Hence  from  (4.17)  we  have 

Z  =  lim  VT  [Z<,(i)cos0  +  Zo(i)sin6] 

i  —foa 

r  - 

=  lim  —  t2,(r)cos0  +  Z„(.y)sine] 
j  » 0  s 

=  lim  —  (Z,fs)x  +  Zois)y  ] 

■f  ~-t  0  s 

=  flii  +  biy,  (4.21) 


where  from  (4.18)  and  (4.19)  it  is  easy  to  show  that  the  asymptotic  behaviour  as  r  — >  0  is 

Z,  =aiS  +  0(s^X  Z„  =  bis  +  Ois^). 

Here,  ai  and  (;[  are  to  be  found  numerically.  Substituting  (4.22)  into  (4.13)  yields 


Z  =ki  +k2[axx  +  (1  +  bi)y  ]. 


(4.22) 


(4.23) 


Comparing  (4.23)  to  (4.11),  the  effect  of  the  vortex  is  seen  to  establish  a  new  equilibrium  mixture  mass 
fraction  profile  within  the  core  region. 

Another  quantity  of  interest  is  the  magnitude  of  the  gradient  of  the  mixture  mass  fraction,  defined  as 


(4.24) 


Note  that  G  is  a  function  of  /? ,  <!>  and  Sc .  Substituting  (4.23)  into  (4.24),  we  sec  that  the  magnitude  of  the 
gradient  in  the  center  of  the  core  is  given  by 

G  I  +(1 +/’■)']'  '•  (4.25) 

Figure  11  is  a  plot  of  G  versus  R  for  Sc  =  0  -  ^2  -  I-  Note  that  as  /?  — >  “o,  G  0  which  implies  that 
there  are  no  gradients  within  the  center  of  the  core,  and  hence  no  diffusion  across  it.  Thus,  the  core  is  com¬ 
pletely  established.  From  this  figure  it  is  clear  that  the  core  is  95%  established  m.  R  ~  50.  This  is  in  agree¬ 
ment  with  Figure  10  where  a  well  defined  core  is  clearly  visible  at  /?  =  100. 
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4.3.  LARGE  VORTEX  REYNOLDS  NUMBER  ASYMPl'OTICS 

To  invcstigaie  ihc  local  structure  of  the  model  problem  described  in  Section  4.1  at  large  vortex  Rey¬ 
nolds  numbers  with  R  >>  1,  we  begin  with  the  equation 

Z;  +  4r  11  -  exp(-y?  r^/4n]ZQ  =  (4.26) 

oc  A 

which  is  (4.12)  with  r  scaled  as  r  =  r  R  .  The  above  equation  is  solved  subject  to  the  initial  and  boun¬ 
dary  conditions 

2  =  ki  +  k^y  at  /  >  0,  F  ->  oo,  0  <  6  <  2n,  and  i  =0,  F>0,  0<9<  2n.  (4.27) 

OUTER  INVISCID  SOLUTION  -  Z^ 

In  the  limit  R  ->  «>,  an  outer  inviscid  region  exists  and  is  governed  by  the  following  inviscid  version  of 
(4.26), 

Z^^+^Ze^^O,  (4.28) 

r 

where  Z^  denotes  the  mixture  mass  fraction  in  the  outer  inviscid  region.  This  equation  is  valid  for 
F  =  0(1)  and  t  «  0(R),  so  that  the  effect  of  a  potential  line  vortex  or  an  Oseen  vortex  is  the  same.  The 
solution  is  easily  found  by  the  method  of  characteristics  and  is  given  by 

Z“  =  + /:2  Fsine„,  (4.29) 

where  6u  is  the  characteristic 

(4.30) 

r" 

INNER  SOLUTION  -  Z' 

As  F  ->  0,  ilie  solution  (4.29)  becomes  singular  so  the  viscous  tenns  in  (4.26)  must  be  retained.  We  intro¬ 
duce  the  similarity  type  variable 
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where  6  <<  1  and  will  be  chosen  in  the  course  of  the  analysis.  Defining  to  be  the  inner  variable,  and 
expanding 

2' =  Z{s,Q).  (4.32) 

we  see  that  Z  satisfies  the  convection-diffusion  equation 

‘/^(Z  -5ZJ+  [1  -exp(-/?552/4)]Ze=  V^Z.  (4.33) 

os  oSc R 

To  facilitate  matching  with  the  outer  inviscid  solution  (4.29),  we  change  to  a  new  (non-orthogonal)  coordi¬ 
nate  system  in  which  the  characteristic  Gq  is  one  of  the  coordinates.  Thus,  the  above  equation  in  terms  of  s 
and  00  becomes 

‘^(Z  -jZ,)-  -ijexp(-y?5i^/4)Ze„ 


1  ±  _2_J_  1  A  _2__L 

bSc  R  1  3^  bs^  300  s  Bs  5^^  30o  s^  30o  f 


To  balance  both  sides  we  must  choose 


or,  in  terms  of  5, 


bR  Sc  6^’ 


b=:{R  Scr'‘\ 


For  this  scaling,  exp(-/?  5.v^/4)  is  exponentially  small  in  the  limit  R  and  so  again  the  effect  of  a 

potential  line  vortex  or  an  Osccn  vortex  is  the  same.  The  leading  order  equation  in  a  8-expansion  is  given 
by 

■/2(Z  -xZ,)=  ^  Ze^Oo-  (4.37) 

This  equation  can  be  solved  by  separation  of  variables  and,  after  matching  with  the  outer  solution  (4.29),  the 
solution  is  found  to  be 


Z  =  s  exp  — -  sin  Gq. 
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The  inner  viscous  solution  is  now  given  by 

Z'  =  ki  +  k2  s  exp  — ^  sin 00-  (4.39) 

35” 


COMPOSITE  SOLUTION 

The  composite  solution 

*  ^  ^  ' 

Z=ki+k2rexp  — —  sin  0-:^  ,  (4.40) 

3P  r 

.  J  ^  y 

is  valid  for  times  t  <<5"'.  Figure  12  displays  the  composite  solution  Z  versus  s  for  0  =  ±jt/2,  with 
k  \  =  0.5,  ki-  1,5  =  0.05,  and  t  =  1.  Note  the  agreement  with  the  numerical  solution  given  in  Figure  10. 

The  magnitude  of  the  gradient  of  the  mixture  mass  fraction  can  be  obtained  by  substituting  (4.40)  into 
(4.24).  Then  as  5  0,  the  leading  order  result  yields 

G  =  — -y-  exp  ^  |cos(0  -  .  (4.41) 

5s^  [  3s  6s 

Figure  13  is  a  plot  of  G  versus  s  for  0  =  jr/2,  /t2  =  L  and  5  =  1/10.  Note  that  the  maximum  amplification 
is  about  5"'.  Tlie  limit  /?  ->  0  corresponds  to  G  0,  and  so  the  large  spikes  have  been  created  by  the  vor¬ 
tex. 

The  core  radius  growth  can  be  determined  from  (4.31)  and  (4.36),  and  is  found  to  be 

^ - .  =  constant ,  (4.42) 

or,  in  terms  of  dimensional  quantities  (denoted  by  stars), 

—  =  constant.  (4.43) 

This  is  the  main  result  of  Marble  (1985),  who  showed  constant  =  0.5092  +  0(Vd  /F  ).  Finally,  we  note 

from  (4.43)  that  the  spreading  of  the  viscous  core  due  to  the  vortex  field  is  (F*  /D*  )^^  times  as  large  as 

that  obtained  from  diffusion  alone  in  the  absence  of  the  vortex,  as  previously  pointed  out  by  Marble  (1985) 
and  Karagozian  and  Marble  (1986). 


-15- 


5.  NUMERICAL  SOLUTION 


In  this  section  wc  present  selected  numerical  results  to  the  full  system  (2.1)  subject  to  (2.2).  This  sys¬ 
tem  was  solved  by  a  standard  finite  difference  scheme,  as  previously  describeri  in  Section  3.  To  illustrate  the 
numerical  solution  of  (2. 1-2.2)  of  the  continuous  evolution  from  nearly  frozen  flow  to  near  equilibrium  flow, 
we  produce  here  one  result  corresponding  to  the  case  Sc  =  (])  =  1,  /?  =  28,  Zc  =  30.  Da  =  30,  P  =  1, 

and  Pr  -  1  -  2/30.  This  case  corresponds  to  the  asymptotic  results  of  Section  3  as  displayed  in  Figures  3- 
6.  Figure  14a  is  a  plot  of  the  time  slices  of  the  temperature  profile  7’,  while  Figure  14b  is  a  plot  of  the  time 
slices  of  the  reaction  rate  term  Q.  Time  increases  from  the  bottom  lelt  comer  to  the  top  right  comer.  In 
Figure  14a,  a  hot  spot  develops  initially  within  the  viscous  core  (frame  1)  and  rapidly  develops  into  an  iso¬ 
lated,  almost  circular  flame  which  grows  as  time  increases  (frames  2-8).  At  a  later  time  (frame  7)  two 
diffusion  flames  are  clearly  visible  at  the  edges  of  the  plot,  and  as  time  increases  further  (frame  8)  they 
move  towards  the  viscous  core.  Eventually  the  two  diffusion  flames  will  merge  with  the  expanding,  almost 
circular  flame  located  in  tlie  viscous  core  region.  Figure  14b  shows  corresponding  results  for  the  reaction 
rate  term.  Initially  (frame  1)  a  single  isolated  point  in  the  rate  term  is  visible.  As  time  increases  (frames  2- 
8)  this  single  point  develops  into  an  almost  circuha  ring  and  spreads  in  lime,  leaving  within  it  the  burnt 
core.  The  reaction  rates  of  the  two  diffusion  flames  arc  given  in  frames  (7-8). 


6.  CONCLUSIONS 


The  distortion  of  flames  in  flow's  with  vortical  motion  has  been  examined  by  means  of  asymptotic 
analysis  and  numerical  simulation.  The  model  consisted  of  a  constant-density,  one-step,  irreversible 
Arrhenius  reaction  between  initially  unmixed  siiecies  occupying  adjacent  half-planes  which  were  then 
allowed  to  mix  and  react  in  the  presence  of  a  vortex.  The  coiiunuous  evolution  of  itie  temperature  and  mass 
fraction  fields  from  initially  unmixed  to  near  equilibrium  flow  was  Ibllowed.  Emphasis  w'as  placed  on  the 
ignition  time  and  location  as  a  function  of  vortex  Reynolds  numlier  and  iniual  temperature  differences  of 
the  reacting  species.  In  the  ignition  regime,  the  case  of  near  equal  initial  temperatures  was  considered.  The 
effect  of  increasing  the  initial  temperature  ratio  P/  slightly  from  unity  was  to  enhance  ignition  for  any  fixed 
vortex  Reynolds  number  R ,  while  decreasing  P,  slighUy  from  unity  had  the  opposite  effect.  However,  for 
fixed  P/-,  increasing  the  vortex  Reynolds  number  from  zero  had  little  effect  on  the  overall  ignition  time.  The 
ignition  location  occurred  in  the  region  of  the  initially  hotter  reactant  for  small  R ,  and  was  seen  to  spiral 
clockwise  towards  the  viscous  core  center  as  R  increases.  Finally,  numerical  solutions  of  the  ignition  equa¬ 
tions  indicated  that  the  hot  spot  was  both  axisymmeiric  and  diffusion  controlled  for  any  non-zero  vortex 
Reynolds  number.  Thus,  the  asymptotic  analysis  of  Dold  (1985,  1989)  could  be  used  in  describing  this  l(x;al 
structure.  The  case  of  different  initial  temperatures,  and  the  future  time  development  of  the  hot  spot,  will  be 
considered  in  a  future  manuscript. 

In  the  diffusion  flame  regime,  a  simple  mixlel  problem  was  proposed  to  investigate  the  effect  of  the 
vortex  on  the  flame  sheet.  The  model  consisted  of  a,ssuming  that  tiie  time  scale  of  the  vortex  was  faster  tlian 
that  of  diffusion  alone,  so  that  on  this  time  scale  appropriate  asymptotic  expansions  were  used  in  describing 
the  interaction  about  the  flame  sheet.  This  model  problem  allowed  for  the  full  vortex  flow  field  without 
making  the  strain-shear  approximation  of  Marble  (1985).  Many  of  the  conclusions  of  Marble  (1985)  were 
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vcrified,  including  the  similarity  rule  for  the  core  radius  growth  in  the  limit  of  large  vortex  Reynolds 
number. 
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Table  11.  IGNITION  TIMES  vs  VORTEX  REYNOLDS  NUMBER  (R) 
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Ti  at  IGNITION 


Figure  3.  Surface  plol  of  (a)  Ti  and  (b)  Fi_i  for  R  =  28.0,  (ii  =  -2,  and  Sc 


INERT  SPECIE  (  Fi,i  )  AT  IGNITION 
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higurc  7.  Plot  ol  Uie  igiiiiioa  lucaiioas  yig  versus  for  various  values  of  Ute  voriex  Reyaulils  number 
wiili  Sc  =  (|»  -  1.  -  -2.  uial  using  ilie  poieniial  voriex  dcliacd  by  V  =  RSc/r. 


figure  8.  Scheiuuiic  sliuwiiig  ihc  assumed  insumtaacous  flow  Held. 
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Figure  11.  Flui  of  G  versus  vorlex  KcynolUs  Huiiibct  for  Sc  -  ^  -  k^  -  1. 


TIME  SLICES  of  REACTION  PRODUCT 
(BETAx=-2,  R=28) 
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Figure  14.  Plot  ol  the  time  slices  of  (a)  Uie  leinpcrulure  profile  T  arid  (,b)  die  reaction  rate  tenn  il.  Time 
increases  from  die  bottom  left  corner  to  top  right  corner. 


TIME  SLICES  of  TEMPERATURE 
(BETA-p=-2,  R=28) 
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i  igure  14.  Plot  of  tlie  time  slices  of  (a)  ifte  temperature  profile  T  and  (b)  die  reaction  rate  tenn  U.  Time 
ncrcascs  from  the  bottom  left  corner  to  top  right  corner. 
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